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ABSTRACT 



A variable conductance heat pipe was designed and con- 
structed to study the performance characteristics of the 
gas loaded heat pipe with particular emphasis on the effect 
of gravity. Heat inputs were varied from 10 to 50 watts 
and the operating position of the heat pipe was varied from 
horizontal to vertical with the condenser above the evapora- 
tor. Condenser temperature profiles are presented for 
various heat inputs and positions for the conventional and 
gas loaded heat pipe. The possibility of stratification in 
a horizontal heat pipe operating with a non-condensible 
gas that has a larger molecular weight than the working 
fluid is discussed based on observations presented. 

Additionally, a one dimensional analysis was formulated 
and used to compare with the experimental results. The 
results show a good agreement between theory and experiment, 
indicating that the developed mathematical model adequately 
describes the physical situation. 
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NOMENCLATURE 



Cross sectional area of vapor space 


ft 2 


Cross sectional area of wick-wall 
combination 


ft 2 


Outside surface area of condenser 
section 


ft 2 


Diameter of vapor space (equals 2R v ) 


ft 


Gravity 


ft/sec 2 


Heat transfer coefficient 


BTU/hr °F 
ft 2 


Latent heat of vaporisation 
Permeability 


BTU/lbm 


Thermal conductivity 


BTU/hr F ft 


Length of condenser section 


ft 


Condensation rate 


2 

lbm/ft sec 


Pressure 


lb/ft 2 


Heat input 


BTU/ft 2 sec 


Electrical resistance 


ohms 


Inside radius of pipe 


ft 


Local vapor space temperature 


°R 


Ambient temperature 


°R 


Voltage 

Mass ratio of working vapor 
to total density 


volt 


Surface tension 
Void fraction of wick 


lb/ft 


Kinematic viscosity 


ft 2 /sec 



p Density 



lbm/ft 3 



y Dynamic viscosity 



lb-sec/ft ^ 



Lennard-Jones force constant 

ftp Collision integral based on Lennard- 
Jones potential function 



ft 



0 K 3/2 



sec 



lbf ft' 



6 Distance between wires in screen 

mesh wick 



in 



SUBSCRIPTS 



a Working fluid 

b Non-condensible gas 

s Stainless steel 

L Liquid working fluid 
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I . INTRODUCTION 



A. BACKGROUND 

The heat pipe is a device capable of transferring large 
amounts of heat while operating at a nearly constant temper- 
ature. The major development of this device has occurred 
within the last decade although the original concept was 
proposed in 1944. This proposal, by R. S. Gaugler [Ref. 1] 
of the General Motors Corporation, was not acted upon at 
that time. A similar proposal by L. Trefethen [Ref. 2] was 
made to the General Electric Company in 1962 but was not 
investigated. The advent of the space program led to a 
renewed interest in the heat pipe in 1964 by G. M. Grover 
et al [Ref. 3] of the Los Alamos Scientific Laboratory. 

The possibility of developing a reliable, efficient and 
lightweight heat transfer device for application in space 
equipment was too appealing to allow the concept to again 
be overlooked. 

B. HEAT PIPE OPERATION 

The sealed heat pipe is lined with a capillary wick 
structure. The vessel is evacuated and partially filled 
with a fluid referred to as the working fluid. During 
operation the working fluid is continuously evaporating in 
one section of the pipe and condensing throughout another 
section (see Figure 1). The condensate is returned to 
the evaporator section by means of capillary force acting 
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through the wick structure. This change of phase process 
results in a nearly constant temperature throughout the 
heat pipe. The near isothermal operation of the heat pipe 
was reported in 19^5 by K.F. Bainton [Ref. ^]. The iso- 
thermal characteristic of the heat pipe allows for a more 
efficient heat transfer device than that achieved with an 
equivalent solid metal rod. An effective thermal conduc- 
tivity for a lithium heat pipe of 1,000 to 10,000 times that 
of an equivalent metal rod was reported by workers at RCA 
in 1967 [Ref. 53- 




Heat pipe vapor flow diagram • 
figure 1 



The introduction of a non-condensable gas along with 
the working fluid alters the performance of the heat pipe 
appreciably. During operation the non-condensable gas is 
driven to the condenser along with the working vapor. 



10 



Unlike the vapor however the gas does not condense and 
thereby forms a high concentration of the non-condensable 
gas in a portion of the condenser. As the concentration of 
the non-condensable gas increases the amount of working 
fluid present and the condensation rate decrease. The 
temperature of the pipe in this portion of the condenser 
will therefore become less than in a portion where the non- 
condensable gas is absent or is a very small amount of the 
concentration. Since the pressure of a saturated vapor is 
directly proportional to its temperature, the amount of 
working vapor present will decrease with the temperature. 
For a large enough concentration of non-condensable gas 
in a portion of the condenser, ambient temperature will 
be approached. The small amount of working fluid present 
forming the partial pressure of methanol will not condense, 
thereby reducing the heat transfer area and the conductance 
of the heat pipe. 

If the evaporator temperature increases, the partial 
pressure of the working fluid increases thereby compressing 
the non-condensable gas into a smaller volume which in 
turn increases the heat transfer area. The pipe's conduc- 
tance is thereby increased. This characteristic of the 
gas loaded heat pipe enables it to have the potential to 
experience smaller operating temperature variations than 
a conventional heat pipe for similar changes in heat input. 
Therefore, the process by which the non-condensable gas 
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Temperature 



collects in the condenser has a direct bearing on the 
heat pipe's performance. 

Two analyses have been developed to give insight into 
the behavior of the gas loaded heat pipe. The first, the 
flat front theory, was developed by W. Bienert [Ref. 6]. 
This theory assumes that during steady state operation of 
the gas loaded heat pipe the vapor and gas are separated 
in the condenser by a vertical flat front. The axial heat 
conduction of the pipe is considered negligible in this 
theory. The possibility of a pressure grac'.ient being 
present in the pipe is excluded and diffusion is not con- 
sidered. With all of these assumptions an analysis of the 
heat pipe performance is straight forward amd proves ade- 
quate in the evaluation of several heat pip-e parameters. 
Specifically the volume of the gas reservoir and the ratio 
of condenser to reservoir volume can be determined. Figure 
2 represents the flat front theory. 



h- 



evaporator 



condenser 




— • — >- 

reservoir 



Vapor-gas flat front model 
figure 2 
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More recently, Marcus and Edwards [Ref. 7] have developed 
the diffuse front theory. In this analysis the contribution 
of binary diffusion between the non-condensable gas and the 
working vapor is considered along with axial conduction 
in the pipe wall. In addition to heat transfer from the 
condenser by convection the contribution due to radiation 
is included. The inclusion of these phenomenon results 
in the vapor-gas interface spreading out into a more 
realistic profile. This analysis develops accurate 
temperature and concentration profiles for the condenser 
section . 

This analysis was used to formulate the TRVJ computer 
program which is capable of calculating, for various gas 
loaded heat pipe geometries, among other parameters: 

(a) the wall temperature profile; 

(b) heat and mass transfer; 

(c) gas reservoir size. 

C. OBJECTIVE 

V/hen a container is filled with two vapors of dissimilar 
molecular weight the vapor of larger molecular weight will 
settle to the bottom portion of the container due to gravity. 
In the operating gas loaded heat pipe the non-condensable 
gas tends to collect at the end of the heat pipe away from 
the evaporator, having been convected there along with the 
working vapor. An inherently unstable situation will 
therefore be present in the gas loaded heat pipe if the 
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condenser is above the evaporator and the molecular weight 
of the non-condensable gas exceeds that of the working 
vapor. The tendency of the gravitational force is to 
oppose the convection of the non-condensable gas, thereby 
spreading out the concentration of the gas. This decrease 
in the gradient of the working vapor concentration is 
accompanied by a decrease in the temperature gradient. 
Therefore, it was the objective of this work to investigate 
the operation of a gas loaded heat pipe with particular 
emphasis on the effects of gravity. 

The heat pipe used was designed to operate horizontally, 
vertically with the condenser above the evaporator or any 
angle in between these two positions. This variation of 
the angle of operation allowed a variation in the axial 
component of the gravitational force. The results presented 
herein are the product of both data obtained from a heat 
pipe designed and constructed for this purpose and a one 
dimensional analysis. 
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II. DESCRIPTION OF EXPERIMENT 



A. OVERALL DESIGN 

As mentioned earlier the axial conduction of the heat 
pipe wall tends to spread out the gas vapor interface. To 
reduce this, the heat pipe was constructed from a stainless 
steel pipe with a wall thickness of 0.0635 inches. Fins in 
the form of Acme threads were cut on a lathe along the 
condenser section to a depth of 0.0255 inches. 

One of the two end plates was soldered in place while 
the other was designed to be easily removed from the condenser 
end of the heat pipe. The removable end plate, as shown in 
Figure 3> was in the form of an internally threaded cap. The cap 
was fastened to a threaded flange which was welded to the end 
of the heat pipe. An 0 ring was placed in the flange to 
effect a pressure and vacuum seal. The condenser end cap was 
designed with a 0.25 inch diameter tube soldered to its 
external side. Through a T connection the tube provided the 
necessary inlet for a pressure transducer, vacuum valve and 
means of evacuating and filling the heat pipe. 

A mesh stainless steel screen around the inside wall of 
the pipe was used for the wick. To retain the wick firmly 
against the wall a longitudinal stainless steel spring was 
inserted along the entire length of the pipe. 

The evaporator heating device was a resistance heater 
made of .125 inch wide, flat Nichrome strip powered by a 
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direct current power supply, 
design of the heat pipe. 



Figure 4 shows the overall 
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Overall heat pipe design 



B. DESIGN CONSIDERATIONS 



The performance of a heat pipe is affected by any 
chemical reaction or decomposition of the working fluid or 
the corrosion or erosion of the wick, springy or heat pipe. 

An example is the reaction of water and stainless steel 
resulting in the release of hydrogen gas. In designing a 
heat pipe the materials used for the pipe and the wick must 
be chosen to be compatible with each other and the working 
fluid. For this reason, the screen mesh, pipe and spring 
wire were all made of the same material, stai'.nless steel. 

Because the heat pipe operation is dependent on capillary 
pumping the working fluids surface tension ha.d to be consid- 
ered. In addition, the gravitational force and flow losses 
made the consideration of the working fluid* :s density and 
viscosity important. Methanol, whose properties are slightly 
less desirable than water was chosen because of its good 
wetting characteristic and high compatibility with stainless 
steel. Ammonia was rejected because of its ttoxicity. 

Following Marcus [Ref. 7], the screen mesh size and 
dimensions of the pipe for a given working flluid can be 
related by a pressure force balance equation. This equation 
is based on the principle that the net capiilary head in the 
wick must be greater than the losses due to pressure drop in 
the vapor and liquid plus any losses due to body forces. 

This force balance equation is expressed as “ 
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net body liquid 

capillary - force _> pressure 
head head drop 



vapor 

pressure 

drop 
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and the vapor side pressure drop is assumed to be much less 
than the other pressure components and therefore not 
considered . 

The heat load was anticipated to be in the range of zero 
to 100 watts with the most probable loading between 10 and 50 
watts. The stainless steel screen mesh was available in 
various sizes up to 200 wires per inch. Using this infor- 
mation and the above force balance equation an iterative 
design ensued for an operating temperature of 120 to 170 °F. 
Various screen wire spacings, lengths of pipe and number of' 
layers of mesh were combined with various pipe diameters in 
seeking the proper design. The design centered on providing 
for a pipe from one to two feet long and with a diameter in 
excess of one inch if possible. 
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With a margin of safety, the design w as tentatively set 
for a two-inch diameter pipe of overall length 15 inches with 
5 layers of 150 mesh screen. This length allowed for a 4 
inch heater section and a one inch adiabatic section. This 
was the tentative design because these parameters were for 
a heat pipe not loaded with a non-condensable gas. To 
achieve an optimum design for the gas loaded heat pipe the 
design program developed by Marcus [Ref. 7] was employed. 

In order to use this program the non-condensable gas had 
to be selected. The choice centered around the inert gases 
since they were chemically inactive and provided a good 
spectrum of molecular weights. To assist in emphasizing the 
effect of gravity on the vapor-gas mixture a large difference 
in the molecular weights of the working fluid and the non- 
condensable gas was needed. The choice of krypton allowed 
a molecular weight ratio of krypton to methanol of approxi- 
mately 2.5* 

With the heat pipe parameters tentatively determined, 
Marcus' computer program was used to first determine the 
additional length of pipe needed for the non-condensable gas 
and then to give a performance prediction. Two inches were 
added to the condenser to allow for the proper performance 
of the heat pipe with the non-condensable gas. A check of 
the preceding force balance equation used in the design 
process indicated the heat pipe would still operate satis- 
factorily with the extra 2 inches of length. 
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C. INSTRUMENTATION 



The instrumentation of the heat pipe included both 
temperature and pressure sensors. Seventeen thermocouples 
were placed in line in the grooves of the fins at an interval 
of approximately 0.650 inches on the top of the pipe. An 
additional thermocouple was placed inside the pipe. This 
thermocouple, enclosed in a stainless steel sheath, was 
inserted 2 inches through the evaporator end plate. The 
condenser thermocouples were 0.010 inches in diameter copper 
and constantan wires welded with a bead at the end. The 
bead was made with a Dynatch thermocouple welder. These 
beads were attached to the grooves of the condenser with a 
Unitek model 1065 welder. The other end of all the thermo- 
couple wires were led to a Jones strip and from there to a 
Hewlett Packard 2010C Data Acquisition System. 

Several layers of Johns Manville MinK insulation were 
wrapped around the heater to minimize heat losses. Two 
additional thermocouples were then placed across a layer of 
insulation. By knowing the temperature difference across a 
layer of insulation as well as the thermal conductivity and 
thickness of the insulation the heat loss from the evaporator 
could be estimated. 

To determine the power input to the heater accurately a 
known resistance was placed in series with the Nichrome 
heater strip. By measuring the voltage across the known 
resistor the heater current was determined. This current 
was used with the voltage reading across the heater to 



21 



determine the power input at the evaporator section. A 
Hewlett Packard digital voltmeter was used for these 
measurements . 

The pressure was measured by a Pace KP-15 transducer and 
recorded on a Hewlett Packard 7702B recorder. Figure 5 shows 
a schematic diagram of the heat pipe and its associated 
equipment. Included in this diagram and in Figure 6 is the 
apparatus used to fill the heat pipe with the required 
amount of methanol and non-condensable gas. The filling rig 
used was a previously existing device used by J. Woodard 
[Ref. 8] which was modified to include an inlet for the non- 
condensable gas and an additional vacuum valve. 

D. CALIBRATION AND ERROR ANALYSIS 

The- thermocouples were calibrated in a constant tempera- 
ture bath with a platinum wire resistance thermometer and a 
Rosemont commutating bridge, while the thermocouple readings 
were made on a Leeds and Northrop K-3 precision potentiometer. 
The thermocouple readings were found to vary from the com- 
mutating bridge temperature readings by no more than ±0.11°F. 
The data acquisition system was then checked against the Leeds 
and Northrop potentiometer and found to have an error bound 
of ± 0.4°F. The resistor placed in series with the evaporator 
heater was also calibrated on the Rosemont commutating bridge. 
The resistance was read as 2.013 ± 0.0002 ohms. The Harrison 
pressure recorder and Pace transducer were calibrated with a 
Wallace and Tiernan gauge. The error for the transducer and 
recorder was ± 0.4 psia. 
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An error analysis was performed for all the data 
recorded as well as the other parameters used in computations. 

A summary of the error bounds is shown in Appendix A along 
with sample calculations. The data points on the graphs were 
drawn to include the error bounds of the information presented. 

E. EXPERIMENTAL PROCEDURE 

Prior to proceeding with the experiment, the amount of 
working fluid to be used was determined. The equivalent 
liquid volume of the working vapor for a heat input of 
50 watts was added to the void volume of the wick structure. 

The volume of methanol required was approximately *J5 milli- 
liters depending on the safety factor used. Operation of 
the heat pipe with 30, *10, and finally 55 milliliters of 
methanol produced no discernible difference in the temperature 
and pressure readings. Since wick dry out, caused by a 
shortage of working fluid in the wick structure, can cause 
a breakdown in the capillary pumping process the choice was 
made to use 55 milliliters of methanol. An excess of 
working fluid would not produce any detrimental effects on 
the heat pipe's performance. The krypton gas used in the 
experiment was research grade while the methanol working 
fluid was spectrophotometric quality. 

Prior to assembling the heat pipe, the pipe, wick, end 
cap, and spring were all scrubbed with acetone and rinsed 
three times with distilled water. The assembled heat pipe 
was subjected to the same process prior to attaching the end 
plate. The procedure used in preparing the assembled heat 
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pipe for operation started with evacuation of the heat pipe 
and filling assembly. The filling assembly, in addition 
to supplying the working fluid and non-condensable gas, 
provided for the connection of the heat pipe to the evacua- 
tion pump. A total of approximately 30 inches of 0.25 
inch inside diameter vacuum hose connected the heat pipe 
through the filling assembly to the evacuation pump. The 
pump was able to achieve vacuum, read at the pump, of 
2 x 10" 6 millimeters of mercury. 

When the evacuation was completed the heat pipe and 
filling assembly were isolated from the evacuation pump. 

A measured amount of methanol was then released into the 
heat pipe. The pressure was recorded when the temperature 
inside the pipe equalled ambient temperature . The krypton 
was next released from its high pressure cylinder and allowed 
to flow into the heat pipe, while monitoring the pressure 
and internal temperature. The amount of krypton present 
could then be calculated by using the temperature and the 
difference between the pressure reading and the partial 
pressure of methanol with the perfect gas equation. 

When the internal temperature of the pipe was equal to the 
ambient temperature, it was recorded along with the pressure 
in the heat pipe. 

The heat pipe was then isolated from the filling assembly 
by means of the vacuum valve on the stem of the heat pipe's 
end cap. The pipe was then positioned for the desired angle 
of operation and the heat input was applied. 
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The temperature and pressure readings were monitored 
until a steady state operation was detected. The steady 
state condition was based on the condenser temperatures 
rate of change being less than 0.5 °F per hour. 

The temperatures and pressure were recorded at steady 
state, then the heat pipe was either repositioned or the 
heat input was altered to initiate a new operating mode. 
This process was continued until sufficient data had been 
obtained. 

Shutdown of the pipe required only the removal of the 
heat input and subsequent release of the internal fluids 
when the heat pipe had achieved ambient temperature. 
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III. EXPERIMENTAL RESULTS 



A. CONVENTIONAL HEAT PIPE 

The heat pipe was initially operated as a conventional 
heat pipe with only methanol present. Heat inputs were 
applied from 10 to 50 watts in 10 watt increments in each 
of three heat pipe positions; horizontal, condenser raised 
*J5°, and vertical. The data presented in Figures 7> 8, and 
9 is representative of all the condenser temperature readings 
for the various heat inputs and positions mentioned. As 
these Figures show, the entire condenser operated in an 
isothermal manner as expected. 

Figures 10 and 11 show the evaporator temperature and 
condenser pressure for a particular heat input. The data 
points so obtained for each of the heat pipe positions are 
connected with smooth curves. The evaporator temperatures 
were used to determine the saturated vapor pressure of 
methanol from the data on methanol presented by J.M. Smith 
[Ref. 9]. The pressures obtained from J. M. Smith's table 
were all within ± 5$ of the corresponding pressure readings 
of Figure 10. Therefore, the assumption of operation of 
the heat pipe with a saturated working vapor is confirmed. 

The temperature readings used to enter J. M. Smith's table 
were from the evaporator while the pressure reading used 
for comparison were taken from the end of the condenser. 

The fact that the pressures are equivalent, within the 
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figure 9 
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accuracy of this experiment, and were read at opposite ends 
of the pipe indicates that there is no pressure gradient 
throughout the conventional heat pipe during operation. 

Figure 11 indicates the heat input to the heat pipe 
as a function of the evaporator temperature for the horizon- 
tal and vertical positions. The variation of the evaporator 
temperature for the changes in heat input is shown to be 
linear. The slope for the vertical position in Figure 11 
is less than that for the horizontal position. This indi- 
cates higher evaporator temperatures for the same heat 
input in the vertical position when compared to the horizon- 
tal position. An explanation for this can be found in the 
heat transfer coefficient of the cylinder which is directly 
related to the Nusselt Number by: 
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where the Nusselt Number equals: 
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where L is the significant length and equals: 



the pipe diameter for the horizontal position, 
the overall pipe length for the vertical position. 

Since the overall length of the heat pipe is approximately 
9 times its diameter the heat transfer coefficient for the 
horizontal pipe exceeds that for the vertical case. For the 
same heat input the decreased heat transfer coefficient 
in the vertical position, when compared to the horizontal 
position, requires a greater temperature difference between 
the heat pipe and the atmosphere. This is reflected in 
Figure 11. 

B. GAS LOADED HEAT PIPE 

After evacuation the heat pipe was filled again with 
55 milliliters of methanol and 2.5 x 10 poundmass of 
krypton gas. The heat pipe was operated over a similar 
range of heat inputs and the same three positions. The 
condenser temperature profile for various heat inputs and 
positions are shown in Figures 12, 13 and 14 . 

Figures 12 and 13 representing the horizontal and ^5° 
position, show only a nomial negative temperature slope 
from the evaporator end of the condenser. This temperature 
profile indicated the concentration of working vapor was 
decreasing very little in the direction of decreasing 
temperature. It was anticipated that a larger temperature 
gradient would be present, reflecting a more marked 
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Horizontal position, condenser temperature, methanol and krypton 

figure 12 
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Length Cinches] 

Vertical position, condenser temperature, methanol and krypton 

figure 14 



segregation of the working vapor and the gas. It was further 
anticipated that this front between the vapor and gas 
would spread out over more of the condenser as the angle 
of the heat pipe increased. 

The vertical operating position, as shown in Figure lk 
had a larger negative temperature slope than that found in 
the other two positions. The vertical position appears 
to indicate the absence of a sharp front but rather a 
mixture of the vapor and gas. The concentration of the 
vapor does decrease more markedly with distance from the 
evaporator in this position. This indicates that inspite 
of the gravitational force and the larger molecular weight 
of the non-condensable gas, the gas does tend to accumulate 
in the condenser at the end away from the evaporator. The 
concentration of the non-condensable gas at no point excludes 
the condensation of the w orking fluid in the condenser. 
Therefore, the gravitational force does effect the variable 
conductance property of the gas loaded heat pipe. 

In comparing Figures 7, 8, and 9 with Figures 12, 13 
and 1*J it is apparent that the change in temperature of the 
condenser for a variation in the heat load was less for the 
gas loaded heat pipe. This is also indicated by the larger 
slope found in Figure 11 when compared with Figure 16. This 
is another feature of the gas loaded heat pipe due to its 
variable conductance. As these two figures indicate, the 
operating temperature of the evaporator is higher, for the 
same heat input, for the gas loaded heat pipe. Figures 10 
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and 15 indicate that, as expected, the pressure is greater 
in the gas loaded heat pipe compared to the conventional 
heat pipe for the same heat input. 

The lack of a significant temperature variation in the 
condenser for the horizontal position was investigated 
further after the data had been collected for all heat 
inputs and positions. 

The heat pipe, in the horizontal position was rotated 
around its longitudinal axis in order to allow the thermo- 
couples to take readings along the bottom of the heat pipe. 

In this position the readings were uniform and on the 
average 12 °F lower than those taken on the top of the 
heat pipe. The bottom of the pipe was expected to be cooler 
than the top due to the natural convection. However, the 
temperature slope on the top of the pipe and the absence of 
such a slope on the bottom indicated that the non-condensable 
gas and methanol may have stratified. That is, the krypton 
may have settled on the bottom of the entire length of the 
pipe . 

If the gas vapor had not stratified the nominal tempera- 
ture gradient present on the top of the pipe would have 
been present on the bottom of the pipe. With stratifica- 
tion, the heavier krypton would settle to the bottom portion 
of the horizontal condenser. This would greatly reduce 
the condensation rate if not completely stop it in this 
section. The bottom of the heat pipe being greater than 



ambient temperature can be accounted for by the circumferen- 
tial conduction of the pipe. 

To further investigate this, temperature sensitive 
liquid crystals were applied over a 2 inch wide black paint 
strip which had been applied to the entire length of one 
side of the condenser. The liquid crystal has the charac- 
teristic of changing from clear to a color over a small 
temperature range. The crystals used were R53- These 
indicated three distinct temperature ranges by three 
different colors as indicated below [Ref. 10]: 



Color 


Temperature 


Red 


50.5 ± 0.5 


Green 


52.1 ± 0.5 


Blue 


53.3 ± 0.5 



With the liquid crystals on the pipe the heat input 
was varied until the general temperature range for crystal 
coloring was achieved. Figure 17 shows two photographs of 
the isotherms observed for the indicated power input. In 
the case of a heat pipe the isotherm represents a line of 
constant working vapor density. The fact that the isotherm 
lines run down the side of the pipe in a near horizontal 
manner indicates the working vapor density varies in the 
same manner. Stratification to some degree may therefore 
be present. 
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The effect of natural convection upon the slope of the 
iostherms has not been fully investigated at this time and 
would provide more insight into this phenomenon. The 
implications of the possibility of stratification will be 
discussed in Section V. 



*»5 



IV. ANALYSIS 



A. THEORY 

The analysis developed in this section is for the con- 
denser section of a gas loaded heat pipe. The development 
is through the application of the continuity and momentum 
equations for one dimensional flow, the energy equation for 
the wall, wick, and fin heat transfer as well as Fick's 
law of diffusion and the perfect gas lav;. The development 
resulted in four first order differential equations which 
were solved simultaneously with the use of a Runge Kutta 
routine. The formulation of this analysis is outlined 
below. It should be noted that as Figure 18 indicates, the 
origin of the coordinate system is at the end of the pipe 
away from the evaporator. 
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Starting with the integral form of the momentum equation 



as applied to the section shown in Figure 19 • 




to 



l F x = / pV(V-dfl c ) 



In the following analysis the viscous term has been left 
out of the above equation. This simplifying assumption is 
based on the Reynolds number being very small in the 
operating range of the heat pipe. For a typical heat load 
of 20 vjatts: 



r d = hr = 50 *° 



where. V = r ^ 

h fg A c p 



V = 



68.26 



(481) l (l/6r (1/20) 



V = .04 ft/sec 
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For fully developed laminar flov, T the friction factor f 
can be represented by 64/Rp . This results in a shear force: 

F = SyirV dx 

Considering the entire one foot length of condenser 
the force is in the range of 10**^ pounds or less. Therefore, 
the shear force can be neglected. 

Referring to Figure (19), equation (1) becomes: 

P " *- P + fx dx -* + PS Sin 0 dx = “P v? + Ep v + 57 (p y ) dx ^ 

[V + ^ dx] + m n D dx 



where P is the pressure in the condenser section, V is the 
vapor-gas velocity, V is the velocity of the vapor, p is the 

ci 

vapor-gas density, g is the gravitational acceleration and 
D is the diameter of the vapor space. 

Dividing by A c dx, and taking the limit as dx approaches 
zero gives: 



(2) - || + pg sin 6 = pV ^ + V + m ^ 



dx dx 



a A. 



The continuity equation for the mixture is 



pVA = m it D dx + [pV + -r- (pV) dx] A 

O Cl 



or 



ttD . d / , r \ n 

m a 7 + 37 (pV) = 0 

C 
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(3) 



Applying Pick's First Lav; to species a: 



( 4 ) 



P V 
K a a 



p a v - 




dv; 

dx 



where is the diffusion coefficient for the vapor-gas 
mixture, p is the density of the vapor and v; is p^/p. 

a. ci 

An alternate form of this equation is : 

(5) v(w - 1) = X*' A 



The pressure in the heat pipe condenser, P, is the sum 
of the partial pressure of the vapor, P , and the partial 

3 . 

pressure of the gas, : 

P = P 0 + P. 
a b 

The perfect gas law can be used for the non-condensable 
gas. This is possible because the operating temperature and 
pressure in the heat pipe are well below the critical point 
of the gas. However, during heat pipe operation the working 
fluid vapor is in a saturated state which precludes use of 
the perfect gas equation. Instead the pressure of the 
working fluid vapor can be represented as a function of 
temperature only. Therefore, 

(6) P » f (T ) + p b A T 

‘b 

or, 

(7) P = f (T) + #- p (1 - w)T 
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where R is the universal gas constant, is the molecular 
weight of the gas and f(T) is an equation based on the 
saturated vapor data presented by J.M. Smith [Ref. 9 ]- 

Consider the section of wick-wall shown in Figure 20 for 
an energy balance. 





figure 20 



Assuming the temperature difference between the wick on 
the inside of the pipe and the outside of the pipe is negli- 
gible in the radial direction: 

(8) q + m h« it D dx = q + dq + h(T-T ) u D dx 
I g 00 o 

where, q is the heat being conducted axially by the pipe, 
m is the rate of condensation of the working fluid per unit 
area, h^ is the latent heat of condensation of the working . 
vapor, h is the heat transfer coefficient, T to is the ambient 
temperature and D q is the outside diameter of the heat pipe. 
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From Fourier's conduction la w: 



( 9 ) 



q = -k A 



dT 
p dx 



where A is the cross sectional area of the wick-wall and 
P 

k is the thermal conductivity in' the axial direction of the 
wick-wall combination. The evaluation of the thermal 
conductivity, k, and the heat transfer coefficient are 
developed at a later point in this section. 

Substituting equation (9) into equation (8): 



( 10 ) 



_ h (m m ^ A 

,2 k L A 

dx i 



m h~ rr D 

fg 

k A_ 



The equations necessary for the development of the 
analysis are summarized below: 



( 2 ) 



dP 



dV 



“ dx + P6 sln 6 = pV ~ + V 



dx 



d(pV) ttD V • 

dx + A w m 
c 



(3) i - k (PV) 



<5) 



( 6 ) 



P = f(T) + (1 - w)T 

Ai b 



( 10 ) 



d 2 T 

dx 2 



= - m 



h~ 7i D , A 

* & + !l (t> _ p ) 

k A k v l °° L A_ 



and, 

(ID 



p = w 
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where S is an equation for p a as a function of T, based on 
the data presented by J.M. Smith [Ref. 93* 

As noted in the above equations, the diffusion coeffi- 
cient,*)^ , has been taken as constant. The diffusion 
coefficient equation for low pressure is based on the 
following equation [Ref. 11]: 




O 1 C M, + 

1.858 10“ J T 1 * 5 [ -V~- ] 

“ll 



P °12 fi D 



This indicates that does vary somewhat with temperature. 
However, the above assumption has been made in order to 
simplify the analysis. A mean diffusion coefficient for the 
temperature range involved was used in the program. 

The following non-dimensional variables were used in the 
above set of equations: 



P 








* 



P 



w 



V 



* 




where , 



% 

Po 



and, 



m r = h 



Q 



fg A c 
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m = 



m 

m 



x = r- 



f = 



f(T) 



Substituting these into the above equations: 



(12) p sin 0 - 



P * 

r ° i _ r / 

L p gL J L( 

K o° dx 



H Q 



1TD2p o h fg 



) 2 ^-3 

J gl- 



[pV 52- + V 



+ [( 



dx 

l \Q 



* d( p*V* ) 

S -I 

dx 



ttD 2 p h _ 
o fg 



,2 1 -| r 4 L -i m V 

; gL JL D J D 



03) m* = "i.lpV) 

dx 



w v ■ Tw^TT 



i " pg p o h fgl) 



^ L Q 



dw 

dx 



* * . . P 0 T 0 - * * 



(15) P = f (T) + [- 



I-iP 
b o 



] P T (1 - w) 



, 2 m* hLA * T , 2 _ - , 

06) ^2 - — - < T - T=) - * L Q ■ 

dx * K p o D k A T 

^ P o 



Substituting for the non-dimensional coefficients in the 
preceding equations: p 

o 

Hi = P 0 g L 
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h = r IS i 2 JL 

2 L „2. J gL 



11 D Po h fg 



H - 

N 3 D 



N i, = 



hL 
k A 



N n = 



iJ L 2 Q 

Dk A T 
P o 



N, 



= ~£ > 



i6ql £ 



N~ = 



N 3 N 6 



R p o T o 
M b P o 

4QL 



the equations become 



* . XT dP* X7 r * 1f * dV* , , r * d(p*V*) 

(17) P sin 0 — N, — jr = Np [p V — w + V — ^ — W — L 

dx dx dx 



) 



( 1 8 ) m = 



N 



1 d / * *x 
T (p V ) 



3 dx 



(19) V* = 7 - --r N„N dw 



(^ry 3 6 7* 



(20) P = f (T) + N 7 (l + w)P T 



( 21 ) 



2 * 
d^T 



* # 



dx 



*2 = Nil <T -T„ ) - N 5 m 



5*1 



Investigating the relative magnitude of the coefficient 
of the momentum term in equation (17) with respect to the 
other coefficients for representative values of the operating 
heat pipe, one finds: 

= r HXIOX3.H1 3 n 2 r 1 -i 

2 L ^70 x 3 . 14 ( 1/6 ) ( 1/31 . 1 ) J L (32.2)(3600) J 

N 2 = 2.4xio -5 



A similar investigation of the other terms indicated 

-5 

that W 2 was at least a factor of 10 times smaller than the 
value of the other non-dimensional coefficients. This led to 
the simplifying assumption that the momentum term could be 
excluded from the analysis without loss of accuracy. The 
computation of this coefficient was retained in the final 
program in order to serve as a check on its relative value. 
The momentum equation then reduces to a force balance : 

(22) p *sin 9 - N-, = 0 

dx 



Differentiating equation (20): 



(23) 



* 

dP 
— ?r 
dx 



+ » 7 (I^xs* ♦ 34 T*J < 
dT dx ' dT dx 



1 * # 

N ? ~ S T 

' w^ 
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Subtracting this from equation (22): 



«*> iq ir sln 0 



K ^ + h cs * + 2^ T ,j 



dT dx 



w 



dT 



dT* 

?r 

dx 



1 * * dv; 

- N 7 S T SWg. 

' v; dx 



j . 

Solving equation (2*!) for — -jr and differentiating: 

dx 

^■[^O. + wd-wJO^SV 

dx 7 ax 



dT* n2 



+ [^- G 3 + v: ( 1 -w)G 5 ](^St) 



dx 



+ [jp 0 X + ( 1 - 2w )Gjj ] 

N 7 1 q dx dx 



+ K t [- "ir sin 0 + — 5T2 

k 1 N 7 t dx T 



w . q dT i 
sm 0 — s-J 



dx 



where , 



G, = 



df 



~w 



1 

If 5f 



dT S T 



d 2 f* * * 
a i s T 



G-, = 



dT 



*2 



# * * 

df dS m * df 

— r — sr T - — sr 

dT dT dT 

ITT - 



S* 



(S T ) 



G i, = 



* * 

1 + dS T 



~F 

T 



fT * 
dT S 



Gr- = 



O X X X 

d^S * , S dS 

dT 2 T dT 

-^2 



(<) 2 

dT 
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Substituting equation (18) into equation (21): 



( 26 ) 



p * M 

2-i-j - n 4 (t -t„ ) + 2t(p V ) 

dx ^ h 3 dx 



Differentiating equation (19) and substituting into 
equation ( 26 ) : 



(27) 



2 * 
d^T 



dx 



T M /m * m % , „ M r dS 1 avf 

ir = N,.(T -T ) + Nc-N^L Jr —7 ^r\ 7 

*2 “ 5 6 L dT ’ t w(w - 1) d * 



avf dT 
?r 

dx dx 



* 



2 (w - 1 ) 0 * , dw N 2 , S* d^w n 

t: 6 t T^ r; w(w - 1) 7^2 J 

w (w - 1 ) dx dx 



Substituting equation (25) into equation (27) : 



(28) 0 10 + G n aWy 

d x d x d x 



- G 12 (^r)‘ + 



■ dv; n 2 
1 l 
dx 



dx 



- G, I, sin 0 ^-jr + G, p- sin 0 ^ 3 - 
dx 1J dx ‘ 



where , 



T 10 



NcN fi * * 

C(w - 1) f-± wS G 1 + (w-l)N 5 NgS ] 



11 



N^N, * 2 Nj-N^ * ,, 0 ,x * 

r 5 6 dS , 5 6 c n +hn ( ^ c p 1 

*• + NZ S °1 + NjN 6 S G„J 



W 



dT 



g 12 " n 5 N 6 ~TT 77 s 

w ( w - 1 ) 
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G 13 


N p-N ^ y. 

[-pwS G 3 + N 5 Ng (w - 1 ) s' G 5 3 


(1 

i — I 

O 


1N 5 6 l S 
H 1 N y vc ^ 


g 15 


Vfi s* 

M 1 N ? t «2 



Simplifying equation (25) 



(29) d 2 w 

dx 


o * * 2 * 

_ r d^T „ ,dT / p dT dw 

G ?n — 5T? h — y?> + dp2 — r — ;r 

dx dX dx dx dx 

* 

- G„_ sin 0 + G„h sin 0 

di dx ' ^ dx' 


where , 




G 20 


2 

[jj- G x + v;(l - w)G 1 j] 


G 2 1 = 


, ,2 

G 3 + w ( 1 - w)G 5 ] 


G 22 


[|^ G x + (1 - 2 w)G 4 ] 


°23 


r ^ ^ i 

l n 1 n 7 


G 2k = 


r 1 W -| 

L N 1 N ? ^ 2 J 
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Equations (28) and (29), alone with their associated 
non-dimensional coefficients, form a set of two, second 
order, simultaneous non-linear differential equations. As 
mentioned earlier these equations were solved using the 
IBM 380/67 computer with a Runge Kutta routine. The entire 
program was put in double precision to avoid as much machine 
round off error as possible. 

The value of the heat transfer coefficient for natural 

convection, h, was calculated based on the Nusselt Number, 

N [Ref. 12]. 
nux 



where , 



N 

nux 



m ( n N ) 

gr pr 



25 



where, m = .59 for the vertical position 

and, m = .525 for the horizontal position 



The Nusselt Number for the heat pipe in the horizontal 
position was found to be 11.59, while for the vertical 
position the Nusselt Number equalled 67-84. These values 
were used in the following equation to find the natural 
convection portion of the heat transfer coefficient, h c : 



To this was added the value of the radiation component 
of the heat transfer coefficient, h r , resulting in the 
following final equation: 
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h 



, Jj 1| 
a e (T m -tJ) 

h c + ( T - T ) 

v OO 



Two values for a condenser temperature of 600°R and an 
ambient temperature of 5^0° R are indicative of the relative 
values used: 



horizontal position h = l.Hl 
vertical position h = 1.00 

In determining the axial thermal conductivity of the 
wall and saturated wick combination, k, the conductivity of 
the saturated wick was first determined and then combined 
with the thermal conductivity of the stainless steel wall 
in a parallel manner. The determination of the saturated 
wick thermal conductivity could be viewed as either a 
parallel or series combination of the wire screen and the 
working fluid [Ref. 7]. The two values obtained are: 

k s K L 

v = , LL - ?Qfi 

series <J>k + ( 1 - <j) ) k T 

S ij 

Vrallel = + ^ k L =5.80 

where 4> is the void fraction in the wick, and is equal 
to .387. 

In assuming the value is that of the parallel case and 
combining this with the thermal conductivity of the stainless 
steel wall, the largest possible value of thermal conductivity 
for the wick-wall combination is obtained. This is what the 
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designing of the fins v;as intended to reduce. The axial 
conductivity used was therefore: 



k 



15.2 



BTU 

hr ft 6 F 



B. RESULTS 

The computer program developed to solve the set of 
simultaneous, first order, non-linear differential equations 
is included after Appendix A. Several solutions obtained 
from this program are included in Figures 21 and 22. The 
solutions were obtained for a specific heat input by assuming 
a temperature for the condenser at x equals zero. As a 
first estimate the condenser temperature from the experimen- 
tal phase for the same heat input was used. The computer 
solution indicates in the output the amount of heat the pipe 
can transfer at this temperature divided by the amount of 
heat specified as being the input. The condenser temperature 
was varied until this quotient equalled one. 

The analytic solutions shown in Figures 21 and 22 
indicate a nearly isothermal condenser. The condenser- 
evaporator interface is only from 0.01 to 0.2/2 greater in 
temperature than the other end of the condenser in these 
figures . 
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temperature, methanol and krypton 



V. DISCUSSION 



A. COMPARISON OP ANALYSIS AND EXPERIMENT 

Figures 12 and 14 present condenser temperature profiles 
based on the operation of the gas loaded heat pipe. Figures 
21 and 22 present the condenser temperature profiles based 
on the analysis developed in the preceding section. The 
positions and heat inputs are the same for Figures 12 and 21 
and the same for Figures 14 and 22. 

The analytic solutions all have higher condenser 
temperatures than the data obtained from the operating heat 
pipe for identical heat input and position. The analytic 
condenser temperatures exceed the operational temperature 
in °F by J.2 to 19 percent. The larger difference between 
the analytic and operational condenser temperature is 
found in the vertical operating position. 

The heat transfer coefficient, h, if increased would 
reduce the condenser temperature since for a given heat 
transfer, Q: 



Q = h(T - T J 



The values of h determined earlier were based on natural 
convection over a symmetric cylinder. The addition of fins 
and the removable end plate with its associated flange on 
the pipe served to disturb the flow over the pipe. This 
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disturbance of the natural convection over the heat pipe 
would have more effect in the vertical position where the 
convection flow travels the length of the pipe. Any such 
disturbance would tend to allow greater mixing in the con- 
vective flow . This would allow a greater amount of heat 
to be transferred at the same condenser temperature or the 
original amount of heat transferred at a slightly reduced 
temperature. This is the effect of increasing the value 
of h in the preceding equation. Increasing the value of 
the heat transfer coefficient used in the analytic solution 
will reduce the condenser temperature obtained. This 
appears to therefore be one of the reasons for the differ- 
ence in the analytic and operational condenser temperatures. 

The analytic condenser temperature was nearly isothermal 
as noted earlier. This again differs from the operational 
condenser temperatures obtained. The condenser temperature 
decreases with increasing distance from the evaporator in 
both the analytic and operational solutions. However, the 
analytic solution did not have as large a decrease as the 
operating heat pipe did for the vertical position. 

The difference in the temperature slopes of the analytic 
solution and operating heat pipe are not of such magnitude 
to require the analytic solution to be discarded. In fact, 
both solutions indicate the absence of any front. Therefore, 
the assumptions of negligible viscous and momentum forces 
in the analysis are justified. 
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The analysis is unable to give any insight into the 
possibility of stratification in a horizontal gas loaded 
heat pipe, except to confirm the lack of any vapor gas 
front. The possibility of stratification in an operating 
heat pipe cannot be overlooked. The design of the gas 
reservoir may require alteration to allow the stratified 
gas to flow in and out with heat input and pressure 
changes . 

The heat transfer characteristic of a gas loaded heat 
pipe with stratification will differ from the flat front 
and diffuse models discussed earlier. The stratified heat 
pipe may produce some of the characteristics of the 
diffuse front model with the front turned 90° in the 
pipe . 

B. RECOMMENDATIONS 

Further investigation with the computer program has 
indicated a heat pipe of different dimensions from that 
used in this experiment may yield more interesting results. 
A longer pipe of smaller diameter with a higher heat trans- 
fer coefficient produces a marked vapor-gas front in the 
analytic solution. With such a diffuse front the rotation 
of the pipe through 90° could be attempted to observe the 
gravitational effect on the vapor-gas front. 

In addition, other inert gases could be introduced into 
the existing heat pipe in place of krypton. Lighter gases 
may allow the vapor-gas front to become more apparent while 
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also providing a comparison for the data in this work 
concerning the possibility of stratification. The use of 
liquid crystals could continue to provide insight into 
this phenomenon. 



67 



APPENDIX A 



SAMPLE CALCULATIONS OF ERROR ANALYSIS 



Q 




W Q ■ V 2 + 



-V 2 2 1/2 

( — p W ) 2 ] 

R 2 R 



where V is the voltage, R is the resistance, is the 
error bound in the resistance and Wy is the error bound 
in the voltage. 

Dividing both sides by Q: 



W 



_Q 

Q 



W 



_Q 

Q 



W 



_Q 

Q 



W 



_Q 

Q 



[«vV 2 + 

^ 10 . 8^5 x 
0.0093 

0.93 % 



W r 2 



• 05)' 



1/2 



+ ( 



. 003 >, 2 
2.013' 



1/2 

] 
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Quantity 
T Temperature 
P Pressure 
Q Heat input 
L Length 
D Diameter 

$3 Diffusion coefficient 
h Heat transfer coefficient 
g Gravity 

h^, Latent heat of vaporization 
a Surface tension 
p Density 

v Kinematic viscosity 
V Voltage 
R Resistance 
0 Void fraction of wick 
m Condensation rate 



Error bound 
± 0.5 F 
± 0.4 psia 
± 0.93$ 

± 0.004 feet 
± 0.0004 feet 
± 1.84% 

± 2 . 8 % 

± 0.01 ft/sec 2 
± 0.03 BTU/lbm 
± 2 x 10~ 5 lb/ft 
± 0.01 lbin/ft 3 
' ± 3 x 10~ k ft 2 /hr 

± 0.05 vo its 
± 0.003 ohms 
± 1.0% 

± 2.4% 



69 



COMPUTER PROGRAM 
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DIMENSION A (13 01) 

SIMPSON'S ONE THIRD INTEGRATION 
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